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A way of circumventing the gauge and infrared problems with effective Wilson coefficients is 
shown. Implications of experimentally measured charmless B decays are discussed. 

I. INTRODUCTION 

In past years we have witnessed remarkable progress in the study of exclusive charmless B decays. Experimentally, 

CLEO Q has discovered many new two-body decay modes B — > rj'K^, rj'K^, tt^K^, tt^K^, tt^K^, p^ir^, ujK^ 

2^* ■ and found a possible evidence for B -^ (j)K* . Moreover, CLEO has provided new improved upper limits for many 

Qs , other decay modes. While all the measured channels are penguin dominated, the most recently observed p^Tr~ 

T-H ' mode is dominated by the tree diagram. In the meantime, updates and new results of many B -^ PV decays with 

r^ p — r], rj' , TT, K and V — Lo^cj), p, K* as well as i? — > PP decays will be available soon. With the B factories Babar and 

rj : Belle starting to collect data, many exciting and harvest years in the arena of B physics and CP violation are expected 

to come. Theoretically, many significant improvements and developments have been achieved over past years. For 

example, a next-to-leading order effective Hamiltonian for current-current operators and QCD as well as electroweak 

penguin operators becomes available. The renormalization scheme and scale problems with the factorization approach 

s^ , for matrix elements can be circumvented by employing scale- and scheme-independent effective Wilson coefScients. 

(yQ , Heavy-to-light form factors have been computed using QCD sum rules, lattice QCD and potential models. 

en ■ 

CN ! II. EFFECTIVE WILSON COEFFICIENTS 

o. 
a^■ 

0^ ' Although the hadronic matrix element {0{p)) can be directly calculated in the lattice framework, it is conven- 

r^ • tionally evaluated under the factorization hypothesis so that {0{p)) is factorized into the product of two matrix 

^H, elements of single currents, governed by decay constants and form factors. In spite of its tremendous simplicity, the 

CIh. naive factorization approach encounters two principal difficulties. One of them is that the hadronic matrix element 

Zh ' under factorization is renormalization scale p, independent as the vector or axial-vector current is partially conserved. 

r* • Consequently, the amplitude Ci{p){0) fact is not truly physical as the scale dependence of Wilson coefficients does not 

• ^H , get compensation from the matrix elements. A plausible solution to the aforementioned scale problem is to extract 

rN [ the p dependence from the matrix element (0{p)}, combine it with the /i-dependent Wilson coefficients to form 

d ' /i-independent effective Wilson coefficients. However, it was pointed out recently in ||] that c^^ suffer from the gauge 

and infrared ambiguities since an off-shell external quark momentum, which is usually chosen to regulate the infrared 

divergence occurred in the radiative corrections to the local 4-quark operators, will introduce a gauge dependence. 

Therefore, this solution, though removes the scale and scheme dependence of a physical amplitude in the framework 

of the factorization hypothesis, often introduces the infrared cutoff and gauge dependence. 

It was recently shown in 0] that the above-mentioned problems on gauge dependence and infrared singularity 
associated with the effective Wilson coefficients can be resolved by perturbative QCD (PQCD) factorization theorem. 
In this formalism, partons, i.e., external quarks, are assumed to be on shell, and both ultraviolet and infrared 
divergences in radiative corrections are isolated using the dimensional regularization. Because external quarks are 
on shell, gauge invariance of the decay amplitude is maintained under radiative corrections to all orders. This 
statement is confirmed by an explicit one-loop calculation in |^. The obtained ultraviolet poles are subtracted in 
a renormalization scheme, while the infrared poles are absorbed into universal nonperturbative bound-state wave 
functions. The remaining finite piece is grouped into a hard decay subamplitude. The decay rate is then factorized 
into the convolution of the hard subamplitude with the bound-state wave functions, both of which are well-defined 
and gauge invariant. Explicitly, the effective Wilson coefficient has the expression 



c''^ = c{n)gi{fi)g2{fif) , (1) 

where gi{fJ.) is an evolution factor from the scale /i to mi,, whose anomalous dimension is the same as that of c(/i), and 
52 (m/) describes the evolution from mb to fif {fj,f being a factorization scale arising from the dimensional regularization 
of infrared divergences) , whose anomalous dimension differs from that of c{fi) because of the inclusion of the dynamics 
associated with spectator quarks. Setting /i/ = Mi,, the effective Wilson coefficients obtained from the one-loop vertex 
corrections to the 4-quark operators Oi have the form: 



cf 



where the matrix f gives momentum-independent constant terms which depend on the treatment of 75 . The expression 
of r for AB = 1 transition current-current operators Oi, O2 is given in 0], while the complete result for QCD-penguin 
operators O3, • • ■ , Og and electroweak penguin operators O7, • ■ • , Oio is given in [M. It should be accentuated that, 
contrary to the previous work based on Landau gauge and off-shell regularization [pi, the matrix f given in Eq. (0) 
is gauge invariant. Consequently, the effective Wilson coefficients (0) are not only scheme and scale independent but 
also free of gauge and infrared problems. 

III. EFFECTIVE PARAMETERS AND NONFACTORIZABLE EFFECTS 

It is known that the effective Wilson coefficients appear in the factorizable decay amplitudes in the combinations 
o-2i = c!^ + "/^Cjf.i and 02^-1 = d^^i -f W'^'^i (* = li ' ' ■ i 5). Phenomenologically, the number of colors Nc is often 
treated as a free parameter to model the nonfactorizable contribution to hadronic matrix elements and its value can 
be extracted from the data of two-body nonleptonic decays. As shown in g, nonfactorizable effects in the decay 
amplitudes of B ^ PP, VP can be absorbed into the parameters af^ . This amounts to replacing Nc in a°^ by 
{Nf),. Exphcitly, 

off _ off , ^ ^cff „cff _ off , ^ ^cff (A _^ ... r.\ l-Ys 

[J^c hi [N°")2^-l 

where (l/iV°*^)i = (l/Nc) + Xi with Xi being the nonfactorizable terms which receive main contributions from color- 
octet current operators O. In the absence of final-state interactions, we shall assume that Xi and hence N^^ are real. 
If Xi ^^^ universal (i.e. process independent) in charm or bottom decays, then we have a generalized factorization 
scheme in which the decay amplitude is expressed in terms of factorizable contributions multiplied by the universal 
effective parameters a"^^ . Phenomenological analyses of the two-body decay data of D and B mesons indicate that 
while the generalized factorization hypothesis in general works reasonably well, the effective parameters a1 2 do show 
some variation from channel to channel, especially for the weak decays of charmed mesons. A recent updated analysis 
of i? ^ Dtt data gives [g| 

Nf {B ^ DTT)r^ {1.8 -2.2), X2(-B^ Dtt) -(0.12-0.21). (4) 

It is customary to assume in the literature that {Nc^)i « {Nc*^)2 • ■ • « {N°^)io so that the subscript i can be 
dropped; that is, the nonfactorizable term is usually postulated to behave in the same way in penguin and tree decay 
amplitudes. A closer investigation shows that this is not the case. It has been argued in [^ that nonfactorizable effects 
in the matrix elements of {V ~ A){V + A) operators are a priori different from that of {V — A){V — A) operators. One 
primary reason is that the Fierz transformation of the {V — A){V + A) operators 05.6,7.8 is quite different from that 
of {V — A){V — A) operators Oi,2,3,4 and Og^io. As a result, contrary to the common assertion, Nc*^{LR) induced by 
the {V — A){V + A) operators are theoretically different from Nf^[LL) generated by the {V — A)[V — A) operators 
|Q|. Therefore, we shall assume that 

NfiLL) ^ {Nf)^ « {Nf)^ « (TVf )3 « {Nf)^ « (iVf )g « {Nf)^^ , 

Nf{LR) EE [Nf)^ « {Nf)^ « {Nf)^ « [Nf)^ , (5) 



and N^^{LR) ^ N^^{LL) in general. Since the energy release in the energetic two-body charmless B decays is not 
less than that in i? — > Dtt decays, it is thus expected that |x(two — body rare B decay) | < \x{B — > Dtt)\ and hence 
Nf{LL) « Nf{B -^ Dtt) ~ 2. From the data analysis in the next section, we shall see that Nf^[LL) < 3 and 
Nf{LR) > 3. 

IV. ANALYSIS OF DATA 



We have studied in detail the two-body charmless B decays for Bu.d mesons in 
what follows we show some highlights of the data analysis. 



and for i?, mesons in 



In 



A. Spectator-dominated rare B decays 

Very recently, CLEO has made the first observation of a hadronic h ^ u decay, namely B^ -^ p^ix^ |ll|Jl^ 



preliminary measurement yields: 



B{B^ -^ /9"7r±) = (1.5 ± 0.5 ± 0.4) x 10" 



The 



(6) 



The branching ratios of this mode decreases with N°^{LL) as it involves interference between external and internal 
W^-emission amplitudes. From Fig. 1 it is clear that this class-Ill decay is sensitive to 1/N^^ if N^^{LL) is treated as 
a free parameter, namely, Nf^{LR) — N'f^{LL) = N°^] it has the lowest value of order 1 x 10~^ and then grows with 
1/N^^ . We see from Fig. 1 that 0.35 < l/N^^ < 0.92. Since the tree diagrams make the dominant contributions, we 
then have 



1.1 < Nf{LL) < 2.6 from B"^ -^ p^n^. 



(7) 



Therefore, N^^{LL) is favored to be less than 3, as expected. 




FIG. 1. The branching ratio of B~ — > p°7r~ versus 1/A'^^ . The solid (dotted) curve is calculated using the 
Baure-Stech-Wirbel (light-cone sum rule) model for form factors, while the solid thick lines are the CLEO measurements 
with one sigma errors. 



There are two additional experimental hints that favor the choice N^ {LL) ^ 2. First is the class-Ill decay 
B -^ n Lu. This mode is very similar to p^Tr* as its decay amplitude differs from that of lutt only in the 2(03 -I- 05) 
penguin term which is absent in the former. Since the coefficient (03 + a^) is small and it is further subject to the 
quark-mixing angle suppression, the decay rates of cutt and p^Tr are very similar. Although experimentally only 
the upper limit B{B^ -^ tt^uj) < 2.3 x lO'^ is quoted by CLEO ||l|l, the CLEO measurements B{B^ -> K^uj) = 
(1.5+{5:g ± 0.2) X 10^5 and B{B^ -^ h^cu) = (2.5+J5:? ± 0.3) x 10"'"^ with h ^ n, K indicate that the central value of 



B{B^ 



is about 1 X 10"^. This means 0.4 < l/Nf{LL) < 0.6 (see Fig. 3 of g) or 1.7 < Nf{LL) < 2.5 is 



favored; the prediction for Nf{LL) = 2 is B{B^ -^ cJ7r±) = 0.8 x 10"*^ and 1.1 x lO"'^ in the Bauer-Stech-Wirbel 
(BSW) model ||lj| and the light-cone sum rule (LCSR) analysis [|l^ for heavy-to-light form factors, respectively. The 
second hint comes from the decay B'^ -^ tt+tt^ . A very recent CLEO analysis of i?" -^ tt+tt" presents an improved 
upper limit, B{B° —> tt+tt") < 0.84 x lO"'^ |l§. It also imphes that Nf{LL) is preferred to be smaller d). 



B. B ^ <j>K, 4>K* decays 

The decay amplitudes of the penguin-dominated modes B — » (pK and B —t cpK* are governed by the effective 
coefRcients [03 + 04 + as — jl*^? + 0,9 + '^lo)]- Note that the QCD penguin coefficients a^ and 05 are sensitive to 
N'f^{LL) and N^^{LR), respectively. We see from Figs. 6 and 7 of ji) that the decay rates of i? — > (j)K^*^ increase 
with l/Nf{LR) irrespective of Nf{LL). The new CLEO upper limit B{B^ -> (t)K^) < 0.59 x 10"^ at 90% C.L. 
| p2[ implies that 

with N^^{LL) being fixed at the value of 2. Note that this constraint is subject to the corrections from spacelike 
penguin and Vl^-annihilation contributions. At any rate, it is safe to conclude that N^^{LR) > 3 > N^^{LL). 

CLEO has seen a 3cr evidence for the decay B — > (j)K* . Its branching ratio, the average of (pK*^ and pK*'^ modes, 
is reported to be |]l2| 

B{B -^ (j)K*) = i [B{B^ ^ (t)K*^) + ^(5° ^ 0X*°)] = (l.lloj ± 0.2) x 10"^ (9) 

We find that the branching ratio of i? ^ (f)K* is in general larger (less) than that of i? ^ pK in the LCSR (BSW) 
model. This is because T{B -^ pK*) is very sensitive to the form factor ratio x = Af^ {m^)/Af^ ("T'I)) which 
is equal to 0.875 (1.03) in the LCSR (BSW) model. In particular, B{B -^ pK*) = 0.77 x lO'^ is predicted by the 
LCSR for N°^{LL) = 2 and N°^{LR) = 5, which is in accordance with experiment. It is evident that the data 
of B — ^ (jyK and B -^ pK* can be simultaneously accommodated in the LCSR analysis (see Figs. 6 and 7 of Q). 
Therefore, the non-observation of -B -^ (pK does not necessarily invalidate the factorization hypothesis; it could imply 
that the form- factor ratio A2/A1 is less than unity. Of course, it is also possible that the absence of i? ^ (j>K events 
is a downward fluctuation of the experimental signal. At any rate, in order to clarify this issue and to pin down the 
effective number of colors N^^{LR), measurements of B — ^ <j)K and B -^ (f>K* are urgently needed with sufficient 
accuracy. 

C. B ^ rj'K and rjK decays 

The published CLEO results on the decay B -^ ■q'K: B{B^ -^ -q' K^) = (6.5^1:4 ± 0.9) x IQ-^ are several 
times larger than earlier theoretical predictions |18[-p0[ in the range of (1 — 2) x 10^^. It was pointed out in past two 
years by several authors |||,^ that the decay rate of B ^ -q'K will get enhanced because of the small running strange 
quark mass at the scale m;, and sizable 5'f/(3) breaking in the decay constants /s and /q. 

As shown in 0], if N°^{LL) is treated to be the same as N^^{LR), the branching ratio of {B^ -^ v'K) is 
^ (2.7-4.7) X 10"^ at < 1/N^^ < 0.5 and it becomes (3.5-3.8) x 10"^ when the 77' charm content contribution with 
fZ = —6.3 MeV is taken into account. However, the discrepancy between theory and experiment is greatly improved 
by treating Nf{LL) and Nf{LR) differently Setting Nf{LL) = 2, we find that (see Fig. 2) the decay rates of 
B -^ r]'K are considerably enhanced especially at smaU l/Nf{LR). Specifically, B(-B± -^ r]' K^) at llNfiLR) < 0.2 
is enhanced from (3.6 — 3.8) x 10^^ to (4.6 — 6.1) x 10^'°^ due to three enhancements. First, the 77' charm content 
contribution a2Xc '^ now always contributes in the right direction to the decay rate irrespective of the value of 
N^^{LR). Second, the interference in the spectator amplitudes of B^ — > rj'K^ is constructive. Third, the term 
proportional to 2(03 — a5)Xu '^ -|- (03 -I- 04 — 05)^] ''' ' is enhanced when (iV^'^)3 = (iV^^)4 = 2. Therefore, 
the data of B ^ K-q' provide a strong support for Nf{LL) - 2 and the relation Nf{LR) > Nf{LL). The 
mode B — > rj' K has the largest branching ratio in the two-body charmless B decays due mainly to the constructive 
interference between the penguin contributions arising from the {uu + dd) and ss components of the 77'. By contrast, 
the destructive interference for the 77 production leads to a much small decay rate for B -^ -qK. 




1 / N^" ( L R ) 



FIG. 2. The branching ratio of B± -^ rj'K^ as a function of l/Nf{LR) with Nf{LL) being fixed at the value of 2 and 
rj — 0.370, p — 0.175, ms{mh) = 90 MeV. The calculation is done using the BSW model for form factors. The charm content of 
the rj' with f^, = —6.3 MeV contributes to the solid curve but not to the dotted curve. The anomaly contribution to (r;'|s75s|0) 
is included. For comparison, predictions for 7V° {LL) — N^ {LB) as depicted by the dashed curve with fZ — and dot-dashed 
curve with fZ — —6.3 MeV are also shown. The solid thick lines are the preliminary updated CLEO measurements [[111113] 
with one sigma errors. 



D. B± -^ ujK^ and B± -^ p^K^ decays 



The CLEO observation |||| of a large branching ratio for B^ -^ ujK^, B{B^ -> ujK^) = {l-^tl'l ± 0-2) ^ 10-^ is 
rather difScult to explain at first sight. We showed in [Q that in the absence of FSI, the branching ratio of 5+ -^ ujK^ 
is expected to be of the same order as B{B^ -^ p^K'^) ~ (0.5 — 1.0) x 10~^, while experimentally it is of order 
1.5 X 10~^. We argued that _B+ -^ luK^ receives a sizable final-state rescattering contribution from the intermediate 



states K*-Tr",K*-T],K 



*o^ 



,K p'^,K^p which interfere constructively, whereas the analogous rescattering effect 



on _B+ -^ p'^K^ is very suppressed. However, if the measured branching ratio p^K^ is similar to luK^ , then W- 
annihilation and spacelike penguins may play a prominent role. Likewise, the decay mode B'^ -^ K^ K^ is expected 



to be dominated by inelastic rescattering from p^ p 



intermediate states. 
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